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Abstract. The momentum space zero-range model is used to investigate universal 

properties of three interacting particles confined to two dimensions. The pertinent 

'"O ' equations are first formulated for a system of two identical and one distinct particle 

^ , and the two different two-body subsystems are characterized by two-body energies and 

fj ' masses. The three-body energy in units of one of the two-body energies is a universal 

function of the other two-body energy and the mass ratio. We derive convenient 

fvq ', analytical formulae for calculations of the three-body energy as function of these two 

^ ' independent parameters and exhibit the results as universal curves. In particular, we 

C^^ . show that the three-body system can have any number of stable bound states. When 

the mass ratio of the distinct to identical particles is greater than 0.22 we find that 
fT^ ' at most two stable bound states exist, while for two heavy and one light mass an 

r — , increasing number of bound states is possible. The specific number of stable bound 

C^ ' states depends on the ratio of two-body bound state energies and on the mass ratio 

and we map out an energy-mass phase-diagram of the number of stable bound states. 

Realizable systems of both fermions and bosons are discussed in this framework. 
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1. Introduction 

Quantum mechanics in two-dimensional (2D) systems differs quite markedly from the 
three-dimensional (3D) case in many aspects. A particular example is the centrifugal 
barrier which is zero or positive in 3D, whereas in 2D the s-wave barrier is in fact 
negative. This implies that two particles in 2D are at the threshold of binding even 
when they do not interact, i.e. an infinitesimal attraction will produce a bound state 
[H El |3]. This is seen in the famous Landau criterion which says that potentials with 
negative volume integral will produce a bound state for any value of the strength [1]. 
In fact, even when the volume integral is exactly zero a bound state is still present 
O [HI [TJ m m [To]. These results are in sharp contrast to the 3D arena where it is 
well-known that a finite amount of attraction is required to produce bound states. 

A natural question to ask is how these dimensional differences influence few-body 
states with more than two particles. A direction of research has therefore been to 
derive and establish conditions for the occurrence of properties of such systems. Most 
prominent among the findings is arguably the Efimov effect [llj which has achieved a 
unique position as a mathematical anomaly of three interacting particles in 3D. The so- 
called Efimov states generally correspond to excited states that reside at large distance 
and owe their properties to large two-body scattering lengths. The relative properties, 
however, do not even depend on the scattering length but only on the mass ratios of the 
constituent particles which determine for instance the scaling of energy and size from 
one state to the next. These scalings are completely independent of the short- distance 
details of the inter-particle potentials. 

A number of universal relations have been found and studied in 3D. The classical 
examples are correlations between two observables where a one-dimensional relation 
emerges when the potentials are varied. The Phillips plot is one curve for different 
potentials when the neutron-deuteron scattering length is plotted versus the triton 
binding energy [12], the Coester line of saturation density versus binding energy per 
nucleon for nuclear matter [13], and the Tjon line is the correlation between binding 
energies of three and four nucleons [14j. Other much more abundant 3D systems with 
very little model dependence have been found and baptized halo systems [15] . They 
were at first primarily found as states in nuclei but quickly they were also search for in 
molecules, and in particular as helium dimers and trimers [16j . Universal relations have 
been applied to check for Efimov excited states in nuclear halo systems [IT] and also in 
molecules [18] . 

The study of the Efimov effect has seen a revival in recent years within the context 
of ultracold atomic gases after its initial observation in a gas of Cesium atoms in 2006 
|19j . This was followed by a number of experiments using various atomic species and 
a new subfield dubbed 'Efimov Physics' has since developed (see the recent review in 
Ref. [20]). Some of the latest experiments have demonstrated heteronuclear Efimov 
states [21] and also the presence of some intriguing four-body states [22l [23l [24] tied to 
the Efimov trimers l25l [26] . 
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Common to all the studies and experiments mentioned thus far is that they are 
3D. If one considers instead a 2D setup, then the structure of the three-body states 
is very different. In the case of three equal mass particles the tower of states close to 
threshold is absent and one finds just two bound states in the universal limit where the 
particles reside outside the range of the two-body potential [U EZj. Here we investigate 
a system of three particles where two are assumed identical as function of the two-body 
interaction energies and the mass ratios. As we will demonstrate, one can potentially 
have any number of stable bound states when the parameters of the system are varied. 
We therefore find a much more varied picture of 2D three-body systems in the universal 
limit than has previously been considered. 

Two-dimensional physics is important for a variety of problems such as high- 
temperature superconductors, localization of atoms on surfaces, in semiconducting 
microcavities, and for carbon nanotubes and organic interfaces. One particular 
interesting example is the two-component Fermi gas with attractive interactions in 2D. 
Here the pairing instability is caused by the presence of a two-body bound s-wave state 
[28| [29| [30] . The role of three-body states, however, remains elusive. In the ultracold 
atomic gas community there is great interest in producing quantum degenerate gases 
in low dimensions to investigate some of the basic features of a 2D environment. Early 
experiments succeeded in producing quasi-2D samples of ^^^Cs [HI |32l [33], ^^Na [3^ . 
and ^^Rb [35]. Mixtures of ^°K and '^^Rb have been used to produce 2D gases [Ml [37] , 
and more recently samples with two-component gases of ^Li [381 EH] and ^°K [ID] have 
been studied. Heteronuclear molecules ^°K^^Rb have also been trapped in a layered 
stack of quasi-2D pancakes [H] . 

Quasi-2D system are abundant in ultracold atom physics and we thus expect that 
our predictions can be studied in a variety of these setups. In a real experiment 
the confinement to 2D is typically done using an optical lattice. This introduces a 
transverse energy scale, hujo. Below this scale the physics is effectively 2D while it 
becomes 3D at or above taoQ. In this first study we will assume a strict 2D setup, 
which amounts to assuming that we are always working far below fvjJo. In some of the 
limiting cases considered here this assumption breaks down and in those cases the trap 
must be explicitly included to make predictions for real experiments. However, in the 
weak-coupling limit where the two-body bound state energies of all three subsystems 
is small we expect our results to apply directly. The two-body energy must, however, 
be modified to include the trap scale in the manner outlined in Ref. [12]. In addition, 
confinement-induced resonances occur in low- dimensional systems [13] and have to be 
taken into account when considering the possible few-body bound states. We expect 
that techniques for detection of three-body states and resonances in 3D will also be 
useful in the quasi-2D case, i.e. identification of sharp features in the loss rates [20] or 
radio- frequency association [lH [15] . 

The three-body states we consider here consist of two identical particles and a 
third distinct particle. This can be a system consisting of two identical bosons and a 
third particle, or two fermions with at least two internal degrees of freedom (typically 
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hyperfine projections in the case of ultracold alkali gases). As we consider the s-wave 
channel only in our equations, this means that in the case of two identical fermions, the 
interaction is necessarily zero. While there could be interaction in higher partial waves, 
we neglect these in the present study as they are usually much smaller than the s-wave 
interactions. 

We illustrate the generic system under consideration in figure [1^) as a plane with 
two identical A particles and one B particle. The plane is drawn as a pancake which is 
the typical experimental situation. In the picture we assume that the mass of A is smaller 
than that of B but we will consider both cases below. In figure [TJj) we show another 
interesting setup for which our results are relevant, a multi-layer stack of pancakes with 
particles of both kind in each of the layers. In addition to the in-plane three-body 
states that are possible, a long-range force that acts across the layers can also provide 
binding. Our results are then relevant in the universal limit which applies for ultracold 
polar molecules with small dipole moments. This has been discussed in one- [16] and 
three-dimensional systems very recently [17] , and for the 2D bi- and multi-layers in the 
weak O [71 [HI [9l [10] and strongly coupled cases jH] . This 2D multi-layer with dipolar 
particles has been experimentally realized as mentioned above [H]. 

Occurrence conditions and properties of model independent or universal states, like 
in Efimov or halo physics [l9] , are studied using zero-range models where the properties 
of the systems arise from distances larger than the interaction ranges [Sni [51]. The 
behavior is important as a measure against such asymptotic properties, but clearly also 
directly if these properties can be realized in Nature or in laboratories. Both these 
roles now appear more and more in 2D physics. The purpose of the present paper is 
to provide universal energy relations for three-body systems in 2D. These relations are 
completely general and specific applications only require appropriate two-body input 
parameters. 

The structure of the paper is as follows. The introduction is followed by a sketch 
in section II of the method. Equipped with the appropriate equations the numerical 
investigations follow in section III. The stability of ground and excited states are studied 
in section IV. Along the way we reformulate the integral equation to apply to various 
threshold properties. In section V we briefly sum up and formulate the conclusions. 

2. Formalism and Notation 

We consider a two-dimensional AAB system assuming that the interactions depend only 
on relative distances. This can be a three-body system in a single plane as in figure [1^) 
or three particles in two or three layers when long-range interactions are present as in 
figure [lb). In both situations, the three-particle dynamics effectively happens in a single 
plane. We use zero-range interactions as we are interesting in the for model-independent 
universal limit. This simplifies the formulation of the Faddeev equation for the three- 
body bound state as the contact interaction is separable. In the case of long-range 
interaction across different layers, the use of zero-range interactions assumes that the 



Scaling and Universality in two dimensions 



a) 




Figure 1. Schematic illustration of the generic setup considered, a) A three-body 
system in a plane (depicted more realistically as a pancake with ellipsoidal shape to 
indicate potential trap deformation) with two identical A particles and one B particle 
(top view). The picture indicates that the mass of B is larger than that of A. Two- 
body interaction energies, Eaa and Eab, and coordinate axis are also shown, b) A 
multi-layer setup containing multiple A and B particles (side view). In a system with 
long-range interactions (like dipole-dipole forces) three-body states can be formed from 
particles in the same plane and in different planes. Two kinds of three-body systems 
with particles in adjancent layers are indicated by arrows. Coordinate system is shown 
on the left. 



low-energy properties of the true long-range interaction can be described by an effective 
interaction of zero-range similarly to the van der Waals interaction for neutral atoms 
in 3D. This is true for potentials that go to zero at infinity faster than 1/r^ in 2D [52] , 
which includes the dipolar 1/r^ interaction. 

The 2D Hamiltonian for the three-particle AAB system with a pairwise potential 
is 



H = Ho + Vaa + Vab + Vab 



(1) 



where the two A particles are assumed to be identical bosons. The kinetic energy 
operator is 



Ha 



qi 



+ 



v\ 



q| 



Pl 



qi' 



v\' 



'^f^A,AB '^fJ'AB '^l^B.AA '^^-AA '^fJ'A,AB '^^-AB 

where we use Jacobi relative momenta given in terms of rest frame momenta, kj with 
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2 = a, /3,7, as 

q^ = k^ and P7 = /ia/3 ( — ^ -] (3) 

\ma nip J 

where (a,/3,7) is the cychc permutations of the particles {A^A',B) with masses m^ 
and mn. The reduced masses are 11^8 = "^""^'^ and u^aB = "^^ , ."^^ • The contact 
potential in operator form is 

VafS = Xal3\Xal3){Xal3\, (4) 

where the form factor Xa(5{p-y) = 1, depends only on the relative momentum of the pair. 
The two-body T-matrix for negative energies and zero-range potentials for AA and 
AB subsystems is 

Taj{E) = \xa',)ta',{E){xa^\, (5) 

where, the matrix element of the 2D transition matrices are given by (see e.g. [50j for 
the case of identical particles) 



TA-y{E) 



TTlA'm.y ^ I E 

-4:71 — In ' ' 



-1 



(6) 



niA + m^ \ y EA-y 

where ^ = A oi B and Ea-^ is the energy of the AA and AB two-body bound states. 
We adopt units such that h = 1. The singularity of the two-body scattering equation 
in 2D is subtracted by fixing the pole at the two-body bound state (see e.g. [53]). To 
find the bound AAB system we concentrate on the negative energy region but note 
that for positive energy scattering, the analytic extension can be easily performed in 
(jH]). We note that effective range corrections have recently been discussed for three 
identical bosons [54J and will in general shift the binding energies. Such corrections can 
be included in the current formalism and can presumably be tuned in experiment by 
chosing different Feshbach resonances and trapping frequencies. However, since we are 
mainly interested in the overall structures (such as the number of bound states for given 
masses and interaction ratios) we neglect such corrections in this study. 

The bound state wave function \^ aab) decomposed in terms of the Faddeev 
components is 

\^aab) = \^a) + |*A') + \^b) (7) 

where |^^) = Go{E)Vai3\'^AAB) with the resolvent Go{E) = (E ~ i^o)"^ which is 
nonsingular for bound states. The Faddeev equations are written in terms of the 
transition matrix, which is well defined for the zero range potential and given by ([5]). 
We then have 

\^,) = G,{E)T^p (e - -^^) (|vl/„) + |vl/^)) , (8) 

which simplifies through the separability of the potential to give 

\^,) = Go{E)\xa^)\f,) , (9) 
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where (q-yl/^y) = AQ,^(q^, Xapl^ aab) is the vertex function of the Faddeev component of 
the wave function. 

The coupled set of homogeneous integral equations for the Faddeev components of 
the vertices of the wave function for the bound AAB system follows from ([8]) and ([9]) 
and gives (see e.g. [50] and [51] for three identical bosons) 



JB (q) = 2taa E: 



d'k- 



niB + 2mA 2 

q 

AmAmB 
/A(k) 



X 



E.- 



2mAiTiB 



q- 



q 



(10) 



/a (q) = TAB E: 



d'k- 



niB + 2m^ 
2mA (mB + mA 



-q' 



E. 



qz 



d^k 



niA+mB 1,2 
2mA'mB 

/A(k) 



q 



•J 2mAmB 



:iii 



(q2 + k2) - ^k ■ q_ 

where the AA and AB transition amplitudes are calculated for the energies of the 
corresponding subsystems through ([8]). We note that the coupled set of homogeneous 
equations described above gives the 3D equations for a bound system AAB when the 
transition amplitudes and momentum volume are conveniently substituted by their 3D 
forms [T7] . 

The AAB wave function for the zero-range interaction is a solution of the 2D free 
Schrodinger equation except when the particles overlap. In momentum space, the wave 
function is built from the spectator functions, /a (qA) and /b (qs), which are solutions 
of the coupled equations ( [T0|) and ( [TT]) . The relative Jacobi momentum of particle A 
with respect to the center of mass of AB is given by q^ and for B with respect to AA 
by q^. The AAB wave function becomes 

Ia (q^) + Ia (qA') + Jb (qs) 



'^AAB (qA, Pa) 



E. 



2mA+mB „2 

2mA {mA+rriB ) ^^ 



rriA+mB 2 
2mA'm,B "^ 



(12) 



where the relative Jacobi momenta q^ and q^' re combinations of the pair (q^i, p^) 
with Pa the relative momentum between A and B defined in ([3]). 



3. Universal properties 

We want to find the s-wave three-body binding energy, E^, for a system of one distinct 
B particle and two identical A particles. This yields E^ as function of the two masses, 
{mA,mB), and the two-body binding energies, {Eaa,Eab)- These four quantities are 
the only unknown parameters in the set of equations (JTOl) and ( ITTi) which determine E^. 
Therefore E^ must be a function of these four parameters. By using one of the binding 
energies, Eab, as the unit of energy and m^ as the mass unit, we see from ( ITOl) and ( fTTIl 
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that the scaled three-body energy, e^ can be expressed in terms of only two independent 
dimensionless variables, that is 

-E's 



es 



Eab 
where we used the definitions 

me _ ^3 

tua ' Eab 



^ [ Eaa "isA _ t- / 



m). 



m 



and eAA 



E 



AA 



Eab 



(13) 



(14) 



The universal functions, J-'„, are distinctly different for ground and excited states as 
indicated by the discrete subscript n. We shall predominantly focus on the ground state 
{n = 0) but also extract information about the number of bound states for a given set 
of parameters. When the A particles are identical fermions, the interaction in s- waves is 
zero. The limit caa — ;■ is thus the relevant one for identical fermions in the zero-range 
model employed here. We note that this ignores interactions between the fermions in 
the p-wave channel which are usually always negligible in comparison to the interactions 
in the s-wave channel between non-identical A and B particles. 

Using mass unit tua = 1 and choosing the binding energy I-E^bI = 1, which is 



equivalent to rescaling all momenta by q (k) — )■ a/I-EabIq (k) in the set of coupled 
homogeneous integral equations ( ITOl) and (ITTl) . one obtains 



fsiq) 



IT In 



d^k 



-±2q2 + ,3 



-1 -1 



4m 



^AA 



X 



/-4(k) 



63 + ^Cf + k2 



2m 



q 



(15) 



/A(q) 



477- 



m 



m + 1 
d^k- 



In 



m + 2 
2 (m + 1) 

fB{k) 



es + q 



q^ + es 



d'^k- 



l^lrk^ + k.q 

/A(k) 



£3 + 



m+l 
2m 



M 



k^) + ^k . qj 



X 



(16) 



The coupled integral equations, (TT5|) and (TT6|) . give the scaled three-body energies, €3 
as functions, J^ni^AAyfn), of the two independent parameters. We emphasize that Eaa, 
Eab, and E3 are binding energies while 63 and eAA are ratios of binding energies and 
will be referred to as scaled three- and two-body energy, respectively. An analytical 
solution is not available and we shall instead investigate the functions J-'n by numerical 
means. We first concentrate on J-'q and its dependence on energy and mass. 



3.1. Dependence on two-body energy 

The solutions are firmly established for the completely symmetric case of baa = 1 and 
m = 1, where all masses are equal and all pairs have the same binding energy. Then 
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two, and only two, bound three-body states exist with energies given by [U |21 El EZ] 



.(0) 



16.52 and e 



(1) 



1.27, 



;i7) 



where the superscript denote ground and first excited state, respectively. These results 
have been confirmed by several studies using different numerical techniques [551 EHl EZl 
[58] . The two-body binding energies are Eab = Eaa = E2 = — 4/i^/(/ia^) exp(— 27), 
where a is the scattering length, /i is the two-body reduced mass and 7 = 0.57721 is 
Euler's constant. 



200 



160 - 



120 




'AA 



Figure 2. Scaled three-body binding energy, £3 — J^Q{€AA,'rn), as function of the 
scaled two-body energy (caa) for the ground state with m = 0.1, m = 1.0, and 
m -^ 00. States with energies above the 63 — €aa line are stable. The highest-lying 
curve for large £aa corresponding to infinite mass is the function Ao{£aa) in (|22p . The 
curves for all other values of m are located between the m ^ 1 and m — )■ cxd lines in 
the asymptotic region {eaa ^ 15). 



For the ground state, the scaled three-body energies are shown in figure[2]as function 
of the scaled two-body energy, e^A, for several mass ratios, m. The curve for m = 1 must 
go through the point (1.0, 16.52) corresponding to three identical bosons. The straight 
line, €3 = eAA, is inserted to show that the three-body system has a stable ground state 
with binding energy larger than that of the subsystems. Equivalently, for the binding 
energies we have E^ < E2, where E2 is the smallest of the two-body binding energies 
Eab and Eaa (below we will use the notation €2 = E2/EAB for the corresponding scaled 
quantity). A finite number of stable excited states may also exist with energies above 
the €3 = eAA hne. 

All curves in figure [2] increase rather steeply from baa = where €3 must be finite 
since the binding energy of the two pairs of particles, Eab, is finite. The behavior of 
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Table 1. Values of gi{m) iorm^Q.l, 1.0 and 10.0. 




the curves for small eAA is found to be 

£3 = J'o {^AA,m) = J'o(0,m)(l + QoieAA)), for e^A ^ , (18) 

where J-o(0,m) is a function of mass, and Qo{eAA) is a function of energy with 
Qoi^AA = 0) = and a very large derivative at e^A = 0, hke for example ^0(2;) oc x* 
with s less than 1. The limit eAA -^ applies to two identical fermions. The singular 
behavior of the energy in this weak-coupling limit is also found in many-body studies 
of 2D bosonic systems j59l |60] and in two-component 2D Fermi gases |611 162] . 

The curves in figure |2] all increase almost linearly with eAA for large eAA ■ The 
deviation from linearity is very well approximated by a logarithmic modification factor, 
i.e. 

£3 = ^0 {(^AA, m) ^ eAA In ( — he), for eAA -^ 00, (19) 

where Qi (m) is an increasing function of m as seen from the values given in table [H In 
the extreme limit, we have e^ — > eAA- Here the system is simply that of an A A molecule 
with binding energy Eaa and a single B particle with a negligible contribution to the 
energy. 

As the scaled two-body energy eAA increases, the two strongly interacting particles 
contract into one tightly bound dimer entity for all m. In the limit eAA — t- 00 we are 
thus left with a two-body problem effectively. It is known that for any two-body system 
in 2D we always have a bound state [10], which is the ground state found here. The 
three-body binding energy in the limit where \Eaa\ 3> I-EabI, or eAA ^ 1, is expected to 
be approximately the two-body binding energy between the two identical bosons. The 
other contributions become of much less importance and give a weak logarithmic energy 
dependence. All these features are present in the asymptotic parametrization in flTOl) . 

The lines in figure[2]all have qualitatively the same shape. However, sometimes they 
cross each other. This behavior is better appreciated in figure [3] where a much larger 
variation of eAA is shown. When eAA ^ 1, we see that the three-body energy increases 
with decreasing mass ratio, whereas the three-body energy increases with increasing 
mass ratio for eAA ^ 1- This feature of figures [2] and [3] is also easily seen in figure H] 
below. The limit of eAA — )■ on figure H] agrees well with the calculation of one fermion 
and two bosons with resonant boson-fermion interaction given in [57j . 

3.2. Mass dependence 

The mass dependence of J-q is shown in figure H] for the ground state for several values of 
^AA- When m becomes large, the system consists of one heavy and two light particles. 
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Figure 3. Scaled three-body ground state energy £3 = Jo(eAA, m) for large variations 
of EAA with m — 0.10, m — 1.0, and m = 100. 




m 



Figure 4. The scaled three-body ground state energy, £3 = J^o{eAA,'m), as function 
of mass ratio, m, for caa = 0, caa = 0.1, eAA = 1, and eAA = 10. The dash-dot-dotted 
lines show the asymptotic large-rn values. 



By rewriting ( IT5|) and (TT6|) for arbitrary tAA in the limit of large ni, we obtain mass 



Scaling and Universality in two dimensions 



12 



independent equations, i.e. 



/B(q) 



/A(q) 



TT In 



47rln 




/A(k) 



es + k^-k^-k-q' 



(20) 



d'k 



d'k 



/^(k) 



/A(k) 



k ■ q 



ea + Hq' + k^) 



(21) 



Thus, for large values of m the scaled three-body energy becomes m-independent and 
approaches an energy-dependent constant as seen in figure HI The limit can be expressed 
as 

£3 = -^0 (e^A, m) -> 7Q{tAA, oo), for m -)■ oo, (22) 

where J^o{eAAiOo) is the m — )■ oo line shown in figure O The high-energy behavior is 
related to flT9l) where J-o(e^A5C)o) is obtained from Qi{m = oo). The curve J-o(0, oo) 
approaches the structure of two mutually non-interacting light particles each interacting 
with the same heavy particle. The total energy is therefore a sum of the two light-heavy 
particle energies, that is a total energy of 2Eab, as seen for large m and eAA = in 
figure m This can also be considered an accuracy test of the numerical procedure. 

Decreasing the mass ratio towards zero leads to an increasing scaled three-body 
energy which rises rather quickly when m becomes smaller than about 0.2, as seen in 
figure m In the limit of vanishing m, where the system consists of one light and two 
heavy particles, all three-body energies diverge. This is found from the diverging terms 
in flTSj) and flT6|l . which are compensated by a diverging three-body energy. Furthermore, 
infinitely many bound states simultaneously appear. This behavior can be studied in 
the Born-Oppenheimer approximation, where a screened Coulomb-like energy behavior 
emerges [63] in the limit of vanishing mass ratio m. This can be seen in the next section 
in figure El where we find that 63 oc — for m — )■ when saa = 0, that is 



£3 = J^n (0, m) — )■ Bn — , for m — )■ 0, 



m 



(23) 



where i3„ is a constant depending on whether ground or excited state is considered. 

Through inspection of figures [2] and H] we realize that the minimum possible scaled 
three-body energy is found in the limits caa — >• and m — )■ 00. Correspondingly we 
realize that the maximum scaled three-body energy is found in the limits eAA ^ 1 and 
771 — )■ 00. Actually the three-body energy diverges for small m, so this maximum is 
found only if we exclude the very small m region. 

The crossing behavior in figures [2], |3l and H] can be related to the mass dependence 
of £3. Let us focus on two vertical lines in figure H] for mass ratios m = 0.1 and m = 10. 
We notice that J'o (0, 0.1) > J'o (0, 10) and J'o (10, 0.1) < Tq (10, 10). Therefore there 
must exist a scaled energy, caa, such that J-q {^aa, 0.1) = J-q {^aa, 10), corresponding to 
the point of crossing. The same procedure can be followed for all the possible crossing 
points. 
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Stability of the pure three-body system is determined by the scaled energy es compared 
to the thresholds for binding two of the constituent particles. This means that if the 
scaled energy €3 is larger than all thresholds the three-body system is stable. This 
suggests to measure the excess of es over the largest scaled two-body energy that is eAA 
or 1, when eAA is larger or smaller than 1, respectively. For eAA > 1 the threshold is 
increasing with eAA, and for eAA < 1 the threshold is \Eab\ which in our units is equal 
to 1. 
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Figure 5. The relative tliree-body energy, (es — e2)/e2, of ground and first excited 
states for m = f as function of scaled two-body energy caa- We show the scaled 
three-body energy relative to the possible thresholds, £2, of binding the two two-body 
subsystems. Here £2 is the largest value of €aa and 1. For eaa < 1 the threshold is 1 
and for eaa > 1 the threshold is saa- The peaks result from the different thresholds. 
The solid (black) and dashed (red) lines refer to ground and excited states, respectively. 
The 4, arrow indicates the points given in table [31 while the <— indicates the saa = 
limit. 



4.1. Identical masses 

Consider first ni = 1 where we already know that two three-body bound states exist for 
three identical bosons. We show in figure [5] the stability plot for both ground and first 
excited states. Stabihty for both states corresponds to positive values of €3 > eAA > 1, 
and values larger than 1 for eAA < 1- We observe the known ground state stability 
for all scaled two-body energies, eAA- The peaks are an artifact of plotting relative to 
different thresholds. 
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We see in figure [S] how both these states move with respect to the threshold of 
stability. The variation is so large that we need to use log-log scales. For e^A = 1 we 
find the values in ( IT71) . that is 16.52 and 1.27 for ground and excited state, respectively. 
The ground state remains above the thresholds for all values of eAA- However, the 
excited state approaches the threshold of stability for both large and small e^A, that is 
one excited state is present for m = 1 only when the underlying two-body energies are 
in the interval, 1.00 x 10~^ < baa < 2.36. This is detailed in table |3] for three different 
mass ratios. 
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Figure 6. The relative three-body energies, (es — e2)/£2, for eaa = as function of 
TO. The J, arrows indicate the ttj- values where the states successively emerge (ea = 62). 
When TO, —i> 00 we have £3 — > 2e2. 



For non-identical particles with different interactions and masses, the number of 
stable bound states may be completely different. The conditions for existence of excited 
states can be formulated much more conveniently for different values of m. In general, 
the three-body binding energy of any excited state has to be smaller than the smallest 
threshold for binding any subsystem. In figure |6] we see that for eAA = (directly 
applicable for identical fermions) and m = 1 only the ground state is bound. This 
means that the excited state appears for eAA = when es exceeds the threshold of 1, 
and ceases to exist for e^ = eAA when eAA > 1- These conditions refer to the different 
thresholds, eAA and 1, for eAA > 1 and eAA < 1, respectively. 

We can introduce the threshold condition, eAA = 0, in ( TT5l) and ( TT6|) . The result is 
that only one equation remains, i.e. 



fA{q) 




X 
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Table 2. The mass ratio threshold, m^ , for appearance of stable excited states as the 
mass ratio decrease towards zero for caa = 0. 




d^k- 



/A(k) 



£3 + 



m+1 
2m. 



;q2 + k2) + i-k-q 



(24) 



This means that the two identical particles in the limit of vanishing energy, e^A = 0, 
do not feel each other and the distinct particle only feel the potential separately from 
one or the other. This holds for trivially for identical bosons, and for identical fermions 
since we have neglected higher partial waves in the interaction. 



4-2. Excited states 

A survey of the results are shown in figure |6] as function of m. The ground state is 
always stable and approaches 2 as m — )• oo. The excited states appear one after the 
other as m decreases from m = 1 towards zero where the system consists of two heavy 
and one light particle. From 3D physics we are very familiar with this behavior of a 
denser spectrum for such a system. Ultimately, infinitely many stable bound states exist 
in the limit of m = 0. 

The log-log plot is necessary in figure |6]for the general overview but does not expose 
the precise value of m at the thresholds for the appearance which is defined by a positive 
value in figure [61 These thresholds are instead indicated by arrows and more precise 
values of the mass ratios are given in table El We note in both table and figure how the 
excited states appear closer to each other as m, — )■ 0. 

We now investigate how these results vary as function of e^A and m. We know 
from figure El that for m = 10 or m = 0.1 we have only one or three bound states, 
respectively. An interesting question is for which value of e^A will the number of bound 
states change when we keep m fixed. The results are shown in figure [71 for m = 10 and 
m = 0.1. Again the precise thresholds cannot be seen in this figure, and we present 
therefore the values in table [3l We see that for m = 10 only the ground and first excited 
state are stable, while for m = 0.1 the second excited state is also found to be stable. 

From figures [6l and U\ and from table [31 the question arises of how low m should be 
before the next excited state appears. The condition is that the three-body energy es 
has to exceed both caa and 1. In other words, the thresholds where the state begins 
and ceases to exist are equal. Therefore we search for the m values where solutions to 
f lT5|) and f[T6|) exist for which e^ equals the threshold values of eAA and 1. We conclude 
that both ground and first excited states always exist for any value of m, provided the 
two-body energy caa assumes an appropriate value depending on m, as we now explain. 

A higher number of excited states only exists if the mass ratio m is sufficiently 
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Figure 7. The relative ttiree-body energy, (63 — e2)/e2 of ground and first excited 
states for m = 0.1 and m = 10 as function of tlie scaled two-body energy caa- Here 
£2 is the largest value of ^aa and 1. For €aa < 1 the threshold is 1 and for €aa > 1 
the threshold is eAA- The peaks result from the different thresholds. Solid (black) and 
dashed (red) lines are ground and first excited states, respectively, for m = 10. Dotted 
(green), dash-dotted (blue) and dash-double dotted (cyan) lines are ground, first, and 
second states, respectively, for m — 0.1. The J, arrows mark the thresholds given in 
table [H The <— arrows indicate the eAA — limit. 



Table 3. The energy 
mass ratios m. 



intervals for appearance of stable excited states for different 
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small. For example ii in> 0.22 only ground and first excited states can be present. We 
give in table H] the critical masses for appearance of a higher number of bound states. 
These critical masses are indicated with arrows in figure [8] where we show the m — eAA 
phase diagram for the number of stable excited states. As an example of the use of 
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Table 4. The critical mass ratios, nic, above which only Nc bound states can exist. 
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Figure 8. The m — ^aa contour diagram of regions where stable excited three-body 
states exist. In Region I only the ground state is stable. In Region II the ground and 
the first states are both stable. In Region III three or more states are stable. The solid 
(black) and dashed (red) lines show the limits where the first and second excited states 
are stable. The \. arrows mark the critical masses, {mc), for which the next excited 
state appears, table |4l The t arrows indicate the masses for which the next excited 
state emerges when £2 = 0, table [2] 



figure [8], we can conclude that at most ground and first excited state are present for 
m = 1 and m = 10, whereas no more than three bound states exist for m = 0.1. It is 
worth emphasizing that the critical mass for appearance of the third stable bound state 
is at the kink of the dashed (red) curve in the middle of figure [H The threshold is not 
for eAA = as the almost vertical curve otherwise seems to indicate. 



4-3. Realistic systems 

The results in figure M may seem strange; the number of stable bound states decrease 
when one of the two-body attractions increase. This is similar to the number of Efimov 
states in a 3D system. If we start on the side of the resonance where no two-body 
bound state is present (negative scattering length, a < 0, in a typical zero-range model, 
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see for instance figure 1 in [2D]) then the number of bound states is finite (and can be 
zero). When the attraction to the resonance (|a| = oo) the number of states is infinite. 
On the other side of the resonance (a — )■ 0^) the number of bound states will start 
to decrease again to a finite number. For the Efimov states the explanation is that, 
although all energies decrease with increasing attraction, the three-body states catch up 
with the two-body threshold [15]. The result is that three-body states disappear into 
the two-body continuum, and the number of bound states decreases as a consequence. 
In complete analogy, the energies of the present excited three-body bound states also 
decrease with increasing attraction. However, the two-body thresholds decrease faster, 
and the three-body states become unstable as they merge with the two-body continuum. 

In general, figure [8] shows the number of bound states for any set of parameters 
{eAAi fn)- This pair of numbers can be related to any set of two-body energies and masses 
through the scaling relation in ( 1T3|) and f lT4l) . Thus for any given point in figure El we 
find the true three-body binding energy from figure [2l In the Regions I and II we have 
one and two stable states, respectively. Region III collects parameter intervals where 
more than two stable bound states are present. This region could be subdivided by a 
number of curves similar to the dashed (red) curve. Thresholds for eAA = are shown 
but finite values of eAA extend the mass regions as the boundary between regions II and 
III. We expect that 0.07 < m < 0.22 allow up to three while 0.04 < m < 0.07 allow up 
to four stable bound states. Better values are given in table HI 

Some commonly used alkali atoms for ultracold atomic experiments are ^Li, ^°K, 
^^Rb, and ^^^Cs. These mass ratios range from 0.05 to 22. The pairwise interaction is 
usually tunable by Feshbach resonances. Therefore we have to look for variations of the 
dependence on the two-body binding energies. From figure |8] and table ID we see that if 
m > 0.22 at most two bound states exist. Whether the first excited state is present or 
not depends on eAA- li m < 1 and eAA is larger than about 1, only the ground state 
is stable. If m > 1 and eAA is between about 0.001 and 1 also the first excited state is 
stable. When ni is slightly less than 0.22 both 1, 2, and 3 stable states may be present 
depending on eAA- As m decreases below 0.22, an increasing number of stable states are 
possible when eAA is sufficiently small as indicated in figure El 

The 2D structures are not yet routinely made but a number of experimental 
investigations on identical particles have been reported as discussed in the introduction. 
There we also pointed out some modifications expected from the fact that the 
experiments are only quasi-2D. If we consider a layered system with long-range 
interactions as shown in figure [T]d) then particles can be placed in two or three 
equidistant layers. For example, if we consider three layers and place the identical 
particles in the outer layers and the distinct particle in the central layer, then the scaled 
two-body energy eAA is about 0.25 when all particles have the same mass, see e.g. 
[^ ^^. Precisely two bound states are present in the universal regime for this setup. 
This still holds when the mass ratio is larger than 1. However, when the mass ratio 
decreases, more excited states may be present. 

The physics of highly polarized Fermi gases is also interesting in relation to the 
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present study. The problem of a single impurity interacting with a Fermi sea of 
particles in 2D has generated considerable theoretical interest recently [651 EHl l67] . 
Furthermore, fermionic impurities in Bose-Einstein condensates have recently been 
studied in experimentally in optical lattices (although so far only in 3D) [68]. Theoretical 
studies of these kinds of setups usually focus mostly on the influence of two-body bound 
states. However, we expect that interesting spectra of three-body states can occur in 
these systems and it will be interesting to study how a many-body background such as 
a Fermi sea or a condensate can affect the properties of three-body states. 

5. Summary and conclusions 

In the present paper we investigate the three-body problem in two dimensions. The 
aim is to extract universal properties where any potential with similar (observable) 
constraints is able to describe the model-independent results. This only occurs when 
the properties hinges on large-distance behaviour where the details of the basic two-body 
ingredients are unimportant. Zero-range models are then suitable since all properties 
are determined at distances outside the potential. We therefore employ the established 
method of solving the momentum-space Faddeev equations with zero-range interactions. 

We focus for simplicity on a three-body system with two identical and one distinct 
particle. Three different particles would be essentially as easy to solve but the number 
of free parameters would be doubled and the results much harder to display and digest. 
We leave this generalization for a separate future investigation. We first establish that 
the three-body energy in units of one of the two-body binding energies Eab must be a 
function of only two parameters; the ratio of the other two-body binding energy Eaa to 
Eabi and the mass ratio m = mB/mA- We investigate this two-parameter problem as 
function of the reduced energy and mass parameters. When the identical particles are 
fermions their binding energy is zero in the zero-range model, i.e. Eaa = 0. 

The starting point is three identical bosons where the energies of the two stable 
bound states are well-known. The scaled three-body energy is calculated as function of 
the scaled two-body energy for fixed mass ratio and vice versa. The energy dependence 
is the stronger than the mass dependence and always monotonically increasing, and we 
establish logarithmic dependence of the scaled three-body energy for large two-body 
energy. The mass dependence for fixed scaled two-body energy is in general smaller 
but non-monotonic. It is divergently increasing when the mass ratio approaches zero, 
while an energy-dependent and mass-independent constant is approached for large mass 
ratios. 

A number of excited stable bound states can exist. In general when two equal 
masses are not large compared to the third, ground and first excited states exist and are 
stable. These two lowest states are always possible for all mass ratios but, whenever the 
mass ratio is outside the interval 0.22 < m < 1, the first excited state only occurs for a 
rather narrow band of two-body energies. As the two equal masses become heavier the 
number of stable bound states increase towards infinity. This happens when the mass 
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ratio approaches zero and corresponds to the Born-Oppenheimer hmit. We derive a 
number of threshold values of masses and energies for appearance of these stable excited 
states. Finally, we provide an energy-mass phase-diagram of regions for occurrence of 
specific numbers of excited states and discuss applications for realistic systems. 

In conclusion, we have established the two-dimensional universal energy relations 
for a three-body system of two identical and one distinct particle. These investigations 
are of interest as properties of basic quantum mechanical problems extended to two 
dimensions where the behavior differ qualitatively from that of three dimensions. We 
furthermore expect that our results will become directly relevant in the topical studies 
of two-dimensional systems in cold atomic gases. 
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